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Objectives 

This  program  focused  on  the  development  of  a  synergistic  framework  for  model-based  analysis  - 
model  calibration,  sensitivity  analysis  (SA),  uncertainty  quantification  (UQ),  validation  and  verifi¬ 
cation  (V&V),  and  robust  control  design  — -  for  a  range  of  nonlinear  and  hysteretic  smart  composite 
systems.  A  major  component  of  the  program  focused  on  the  development  of  this  framework  in  the 
context  of  hysteretic  macro-fiber  composite  (MFC)  and  shape  memory  alloy  (SMA)  models.  The  goal 
was  to  develop  a  theoretical  and  numerical  framework  and  initial  infrastructure  to  facilitate  mod¬ 
eling,  design  and  control  for  aerospace,  aeronautic,  industrial,  automotive  and  biomedical  systems 
utilizing  smart  materials. 

Accomplishments 

We  focused  on  the  development  and  implementation  of  a  modeling,  model  calibration,  sensitivity 
analysis  (SA),  and  uncertainty  quantification  (UQ)  framework  for  structures  employing  macro-fiber 
composites  (MFC)  and  shape  memory  alloys  (SMA).  As  illustrated  in  Figure  1,  MFC  are  comprised 
of  lead  zirconate  titanate  (PZT)  fibers  embedded  in  an  epoxy  matrix  and  driven  by  an  interdigitated 
electrode  array.  The  resulting  array  achieves  d 33  motion  which  optimizes  outputs  and  is  sufficiently 
flexible  to  bond  to  curved  surfaces.  The  composite  structures  have  dual  actuator  and  sensor  ca¬ 
pabilities  and  can  achieve  kHz  drive  ranges.  As  detailed  in  [10, 13, 19]  and  depicted  in  Figure  2, 
MFC-composite  structures  are  presently  being  investigated  for  shape  modification  and  flow  control, 
vibration  and  noise  control,  health  monitoring,  and  energy  harvesting.  However,  these  advantages 
come  at  the  cost  of  rate-dependent  hysteresis,  creep  and  constitutive  nonlinearities  that  must  be 
addressed  in  models,  designs  and  control  algorithms. 

It  is  noted  in  [6]  that  shape  memory  alloys  (SMA)  have  the  highest  work  density  of  traditional 
“smart  materials”  but  typically  operate  at  frequencies  ranging  from  1-100  Hz  due  to  the  thermal 
convection  required  relax  the  actuators.  They  are  being  investigated  for  use  as  chevrons  to  improve 
jet  engine  performance  and  decrease  aeroacoustic  noise,  MAV  design,  and  numerous  biomedical 
applications  due  to  the  fact  that  they  are  biocompatible.  However,  like  PZT-based  actuators,  they 
exhibit  rate-dependent  hysteresis  and  constitutive  nonlinearities  at  essentially  all  drive  levels. 


Interdigitated  Electrodes 


Electrode  pzT 


(b) 


Figure  1:  (a)  Macro- fiber  composite  schematic  and  (b)  d 33  actuation  produced  by  voltage  applied  to 
the  patches. 
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Figure  2:  (a)  MFC  for  shape  modification  and  flow  control,  and  (b)  experimental  configuration  used 
for  energy  harvesting. 


Uncertainty  Quantification 

Uncertainty  quantification  is  pursued  in  two  steps:  (i)  determination  of  densities  and  uncertainty 
bounds  for  inputs  consisting  of  model  parameters,  boundary  or  initial  conditions,  and  (ii)  develop¬ 
ment  of  highly  efficient  methods  to  propagate  these  input  uncertainties  through  models  to  quantify 
response  uncertainties  and  errors  bounds.  The  goal  is  to  develop  the  capability  to  provide  model 
predictions  with  quantified  uncertainties. 

We  employed  a  Bayesian  approach  to  construct  parameter  uncertainties  based  on  measured  data 
since  the  consideration  of  parameters  as  random  variables  with  associated  densities  is  a  fundamental 
tenet  of  Bayesian  inference.  Due  to  the  moderate  dimensionality  of  the  parameter  space  in  the  MFC 
and  SMA  models  (e.g.,  18  parameters  and  36  states),  evaluation  of  posterior  distributions  using 
standard  cubature  formulae  was  not  feasible.  Instead,  we  employed  a  delayed  rejection  adaptive 
Metropolis  (DRAM)  algorithm  whose  stationary  distribution  is  the  posterior  density. 

The  chain  and  density  for  a  representative  SMA  parameter  are  plotted  in  Figure  3  along  with 
joint  sample  points  for  two  parameters.  We  first  note  that  the  density  is  non-Gaussian  so  asymptotic 
methods  limiting  to  Gaussian  distributions  will  be  ineffective.  Moreover,  the  parameters  are  highly 
correlated  which  significantly  limits  available  options  for  uncertainty  propagation. 

Stochastic  Galerkin,  collocation  and  discrete  projection  methods  all  rely  on  the  assumption  that 
parameters  are  mutually  independent  or  that  one  has  a  representation  for  the  joint  density  for 
all  parameters.  These  assumptions  are  violated  for  this  model  as  well  as  many  physical  models 
arising  in  general  applications.  To  propagate  uncertainties  and  construct  prediction  intervals,  we 
thus  sample  by  index  from  the  DRAM  chains  used  to  construct  input  densities.  This  yielded  the 
prediction  intervals  shown  in  Figure  4.  The  application  of  Bayesian  model  calibration  and  uncertainty 
quantification  techniques  to  MFC  are  reported  in  [9,10]  whereas  SMA  results  are  reported  in  [1,3,4], 


(a)  (b)  (c) 

Figure  3:  (a)  Chain  and  (b)  density  for  a  representative  SMA  parameter,  (c)  Joint  sample  points 
demonstrating  correlation  between  parameters. 
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Figure  4:  (a)  Calibrated  model  fit  to  SMA  data  and  (b)  data  with  95  %  prediction  intervals. 


Robust  Control  Design 

The  goal  in  this  component  of  the  program  was  to  use  model  uncertainties,  quantified  using  the 
previously  discussed  Bayesian  and  sampling  techniques,  to  construct  robust  control  designs  having 
improved  performance.  We  initially  investigated  this  in  the  context  of  the  sliding  mode  control 
(SMC)  design  depicted  in  Figure  5  and  detailed  in  [1,5,15,16].  The  actuators  were  SMA  tendons 
with  strains  e{t)  generated  by  input  temperatures  T(t).  This  was  modeled  using  the  homogenized 
energy  model  (HEM)  framework  in  [6]  with  the  model  inverse  discussed  in  [14],  The  objective  was 
to  track  a  reference  angle  trajectory  9r(t)  by  heating  or  cooling  the  actuator. 

Since  temperatures  are  not  directly  measured,  we  employed  the  observer 

T0(t)  =  - h[T0(t )  -  Too]  +  7u(t)  +  Hx(t )  +  L[6[t)  -  90(t)\ 

=  f(t,  h,  H)  +  Af(t,  Ah,  AH)  +  (7  +  A7 )u(t)  +  L[9(t)  -  0o(i)] 

where  h  =  h  +  Ah,  H  =  H  +  AH  and  7  =  7  +  A7  are  model  parameters  with  means  h,  H,  7  and 
2a  uncertainties  Ah,  AH,  A7  as  illustrated  in  Figure  3(b).  Furthermore,  x(t)  are  summed  phase 
fractions  from  the  model,  L  is  the  observer  gain,  and  u  is  the  control  input  which  ultimately  consists 
of  the  voltage  applied  to  the  actuator.  The  terms 

f(t,  h,  H)  =  -h[T0(t)  -  T^]  +  Hx{t) 

A f{t,  Ah,  AH)  =  -Ah[T0(t)  -  T^}  +  A Hx(t) 


Inverse  HEM 


u»  l:  - 

Observer 

Figure  5:  (a)  Block  diagram  for  the  sliding  mode  control  design. 
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isolate  the  mean  and  uncertain  components  of  the  observer.  For  the  sliding  surface 


s(t)  =  A  j  [T0(r)  -  Tr(r)}dr  +  [ T0(t )  -  Tr(t )], 
use  of  the  Lyapunov  candidate 


yielded  the  control  u(t) 

Ueq(t)  ~- 

v(t)  = 


V(t)  =  * s2{t ) 

=  ueq(t )  +  v(t)  where 

=  3  \-X[T0(t)  -  Tr(t )]  +  fr(t)  -  f(t ,  h,  H)  -  L[0(t)  -  d0 (t)] 
7  L 

- [A  f(t,  Ah,  AH)  +  A^fueq(t)  +  potanh(s)]  . 


The  experimental  performance  of  this  sliding  mode  control  design  is  compared  with  proportional- 
integral  (PI)  control  in  Figure  6.  The  errors  in  the  square  wave  tracking  performance  illustrate  an 
inherent  physical  limitation  of  SMA  actuators  for  high  frequency  applications  since  convective  cooling 
is  a  relatively  slow  process.  As  detailed  in  [1],  the  sliding  mode  design  generally  provides  superior 
performance  without  the  necessity  of  experimentally  tuning  gains.  This  provides  an  illustration  of 
how  uncertainty  quantification  can  improve  robust  control  design.  We  are  presently  extending  these 
techniques  to  alternative  robust  control  designs  and  other  SMA  and  PZT-based  actuators. 


(a) 


(b) 


(c) 


Figure  6:  Tracking  performance  of  the  sliding  mode  control  (SMC)  and  PI  control  for  a  (a)  0.1  Hz 
sine  wave,  (b)  0.2  Hz  sine  wave  and  (c)  0.1  Hz  square  wave. 


Quantification  of  Model  Discrepancy 

The  previously  discussed  model  calibration  and  prediction  interval  results  were  obtained  using 
statistical  models  of  the  form 

Ti  =  f(U,q)  +£i 

where  f(U,q )  denotes  the  parameter-dependent  model  response  and  measurement  errors  et  are  as¬ 
sumed  to  be  independent  and  identically  distributed  (iid)  and  normal  distributed  with  £i  ~  N{ 0,  a2). 
However,  this  assumption  is  often  violated  in  dynamic  applications.  To  illustrate,  we  model  the  thin 
beam  displacement  data  shown  in  Figure  7  using  the  Euler-Bernoulli  beam  model 


P(x ) 


d2w  dw 
~df2+1~dt 

fX  2 

kpV(t)  / 

J  X\ 


4>dx  + 
4>"dx 


YI(x) 


d2w 

dx2 


+  cl(x) 


d3w 

dx2dt 


<j>"dx 
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which  holds  for  all  test  functions  (j)  £  V  =  {<f>  E  H2(0,  L)  \  <j)( 0)  =  <j>'{ 0)  =  0}.  The  density,  stiffness 
and  damping  relations 

p(x)  =  phb  +  pphpbpXp{x )  ,  YI{x)  =  YI  +  YpIPXp(x) 
cl(x)  =  cl  +  CpIpXp{x) 

reflect  the  differing  geometry  and  material  properties  in  the  region  covered  by  the  patch.  The 
reported  density  2700  kg/m3  for  aluminum  yielded  pb  =  phb  =  0.08775  which  is  fixed  to  ensure  that 
the  remaining  parameters  are  identifiable.  The  parameter  set  is 

q  =  [pb,  7,  YIb ,  YIp ,  clb ,  clp,  kp\ 

where  pp  =  pphpbp,Y Ip  =  YpIp,YIb  =  Y I .  CIp  =  CpIp  and  Clb  =  Cbh ■  The  model  response  is 
the  displacement  y(U,q)  =  w(ti,x;q )  evaluated  at  the  point  x  =  128  mm.  Details  regarding  the 
experimental  setup  are  provided  in  [18]. 

Calibration  over  the  time  interval  [0, 1.0]  yielded  the  fit,  residuals,  and  95%  prediction  intervals 
shown  in  Figure  7.  It  is  observed  that  the  residuals,  and  hence  measurement  errors,  are  clearly  corre¬ 
lated,  and  that  the  95%  prediction  interval  does  not  contain  the  correct  percentage  of  experimental 
measurements.  This  is  due  to  unaccommodated  model  errors,  or  discrepancies,  5(ti)  in  the  statistical 
model 

Tj  =  f{U,q)  +  S(U)  +  £i. 

Damping  and  frequency  analysis  established  that  calibration  on  the  interval  [0.25, 1.25]  can  re¬ 
duce  the  model  discrepancy  as  indicated  by  the  residuals  plotted  in  Figure  8(a).  In  this  case,  the  95% 


Figure  7:  (a)  Model  fit  and  (b)  residuals  for  calibration  on  the  temporal  interval  [0, 1.0].  95%  predic¬ 
tion  intervals  and  experimental  data  in  (c)  the  calibration  regime  and  (d)  extrapolatory  prediction 
interval  [2.0,  2.5]  seconds. 
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Figure  8:  (a)  Residuals  for  calibration  on  the  temporal  interval  [0.25, 1.25].  95%  prediction  intervals 
and  experimental  data  in  (c)  the  calibration  regime  and  (d)  extrapolatory  prediction  interval  [2.0,  2.5] 
seconds. 


prediction  interval  shown  in  Figure  8(c)  contains  the  correct  percentage  of  future  data  points,  thus 
extending  the  validation  regime  for  the  model.  The  development  of  statistical  techniques,  such  as 
periodic  Gaussian  process  representations,  to  quantify  model  discrepancies  for  extrapolatory  predic¬ 
tion  constitutes  a  future  research  direction.  Details  regarding  issues  pertaining  to  the  quantification 
of  model  discrepancy  can  be  found  in  [18]. 
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